Electric Circuit Theory and the 
Operational Calculus' 

By JOHN R. CARSON 

CHAPTER IX 

The Finite Line with Terminal Impedances 

So far in our discussions of wave propagation in lines and wave- 
filters, we have confined attention to the case where the impressed 
voltage is applied directly to the infinitely long line. We have found 
that, by virtue of this restriction, the indicial admittance functions 
of the important types of transmission systems are rather easily 
derived and expressible in terms of well known functions, and the 
essential phenomena of wave propagation clearly exhibited. In 
practice, however, we are concerned with lines of finite length with 
the voltage impressed on the line through a terminal impedance Z x 
and the distant end closed by a second terminal impedance Z 2 . We 
now take up the problem presented by such a system. 

Let K = K(p) denote the characteristic operational impedance of 
the line, and y = y(P) the operational propagation constant of the 
line. We have then 

V-Ae-T+BeT, ^ 

I = ^Ae-^-^Bey\ 

where A and B are so far arbitrary constants. To determine these 
constants we assume an e.m.f. E impressed on the line at x = through 
a terminal impedance Zi and the line closed at x = s by a second termi- 
nal impedance Z 2 . At x = s we have therefore 

Z 2 I=V 
whence from (240) 

%*e-y*A - ^eV'B=Ae-y*+Bey* 

and 

B =-\^ n -c-^A (241) 

1+P2 

where p 2 = Z 2 /A". 

1 Concluded from the issue of January, 1926. 
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At x = we have 

V=E-ZJ 

whence 

(l+PiM + tt-pOa-E, (242) 

where p\=Z\/K\. 

From (241) and (242) we get 

1+P2 z. 



A = 



(H-Pi)(l+P 2 )-(1-Pi)(l-Pa)« 



- n.V-2^ 



and finally 



R _ -q-p 2 ) g - 2 ^ E 

^~(1+Pl)(l+P 2 )-(1-Pl)(l-P2)6- 2 ^ 



j - - fIz- - i r — • (243) 

A+Zi 1 — pi 1 — Pi c _ 2ys 

1+pjl+pa 

If we replace £ by a unit e.m.f. we get the operational formula for 
the indicial admittance A x \ thus 

^\ g-v+nr'*-* _ l (244) 



where 



X=X/(A+Z 1 ), 

_l-pi_ A-Zi 
Ml 1+ Pl X+Z,' 

l-P2 _- £-Z 2 

i+p 2 a~+z 2 " 



P2 = 



K,y,Zi,Zi,ni and M2 are, of course, functions of the operator £. 

The integral equation corresponding to the operational formula 
(244) is 

pmr T °-" A + t)dL (245) 
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Now by (244) we can expand ±1/Z x (.p) ; it is 



Z x {p) K ' '^ K 

e -y(2s+x) e -y(is-x) 
+ X/Xi/X 2 ^ 1"XMiM2 2 ^ 

(24G) 

0-7(4*+*) 

4-X/ZXV2 2 — ^ + .... 

Now we observe that e~ yx /K is simply the operational formula for 
the indicial admittance at point x of an infinitely long line with unit 
e.m.f. impressed directly on the line at x = 0. This will be denoted 
by a x (f). Similarly e~ y( - 2s ~ x) /K is the operational formula for the 
indicial admittance at point (25 — x) with unit e.m.f. impressed directly 
on the line at x = 0. This will be denoted by a-i s - x {t), etc. 

Recognition of this fact allows us to derive a formal solution in 
terms of a series of reflected waves. For let a set of functions 
v ,vi,V2,v 3 , .... satisfy and be defined by the operational equations 

flo = A(£)=X 

V2=\fiifJ, 2 

Va = ^nif4, etc. 
It then follows from the preceding and theorem II that 



(247) 



*-ir 



, , V (t — T)a x (T)+Vl(t-T)a2s- x (T) } 

+v 2 (t—T)a 2s + x (T) + 



(248) 



If, therefore, we know the indicial admittance of the infinitely long 
line with unit e.m.f. directly applied and if we can solve the operational 
equations (247), then A x {t) is given by (248) by integration. This 
solution may well present formidable difficulty in the way of com- 
putation. It is, however, formally straightforward and the numerical 
computation is entirely possible, the only question being as to whether 
the importance of the problem justifies the necessary expenditure of 
time and effort. Without any computations, however, the solution 
(248) admits of considerable instructive interpretation by inspection. 
The first term represents the current at point x of an infinitely long 
line in response to a unit e.m.f. impressed at x = through an im- 
pedance Z\\ v = v {t) is the corresponding voltage across the line termi- 
nals proper. The second term is a reflected wave from the other 
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terminal due to the terminal irregularity which exists there. The 
third term is a reflected wave from the sending end terminal, etc. 
The solution is therefore a wave solution and is expanded in a form 
which corresponds exactly with the sequence of phenomena, which it 
represents. 

The solution takes a particularly instructive form when Z\ = k\K 
and Zi = k*K where ki and k» are numerics. Then 

1 



and 



*-ITEiqRK (249) 

1 i-k 2 1— jfei 
02 i+*i i+k 2 i+fci c " 



I l-ko 

^W-txstI , t , , I- (250) 

. 1 — «1 1 — «2 „ / A , 



i+*i 



If fei =0, & 2 = 1 we have the case of the e.m.f. impressed directly on the 
sending end of the line and the distant end closed through its char- 
acteristic impedance; the solution reduces to 

A x {t)=a x {t) 

as, of course, it should be by definition. 

If £i = and jfe 2 = °° . we have the case of the line open-circuited at 
the distant end, and the solution reduces to 

A x (t) = ]a x (t)-a 2s - x (t)-a 2s+x (t)+a is - x (t) + . . .}. (251) 

Finally, if both k\ and /.'•> are zero, the line is shorted and 

A ,(/) = i «,(/) +a 2s - x (t) +fli f+ *(0 +**-*(*) +••■!■■ ( 262 ) 

The operational equations (247) admit of further interesting and 
instructive physical interpretation without computation. Consider 
a circuit consisting of an impedance Z\ in series with an impedance 
K. Let a unit e.m.f. be applied to this circuit and let v a be the re- 
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sultant voltage across the impedance K. Then, operationally, 

K * 

so that Do, thus defined in physical terms, is the v of equations (247). 
Now let this voltage be impressed on a circuit consisting of an 
impedance 2Z 2 in series with an impedance K — Z 2 so that the total 
impedance is K-\-Z 2 . Let the resultant voltage drop across the im- 
pedance element K — Z 2 be denoted by »ij then operationally 

A' K-Z 2 . 

which agrees with V\ as given by equation (247). 

.Similarly if voltage Vi is applied to a circuit consisting of an im- 
pedance 2Zi in series with an impedance K — Z\ and if v 2 denote the 
voltage drop across impedance K—Zi, then 

V 2 = X/Ul/"9. 

We can thus see physically what the voltages v ,vi,v 2 . . mean in terms 
of simple circuits consisting of K and Z\ in series and K and Z 2 in 
series respectively. 

I shall now work out a specific problem exemplifying the preceding 
theory. The example is made as simple as possible for two reasons. 
First because its simplicity makes it more instructive than when the 
phenomena depicted and the essentials of the mathematical methods 
are obscured by complicated formulas and extensive computations. 
Secondly while the general method of solution illustrated is thor- 
oughly practical we cannot hope to arrive at the numerical solutions 
of the complicated problems without a large amount of laborious 
computations. Problems involving transmission lines with com- 
plicated terminal impedances are among the most difficult, as regards 
actual numerical solution, of any which present themselves in mathe- 
matical physics. On the other hand, the formal solution (248) gives 
at a glance the essential character of the phenomena involved. 

The specific problem we shall deal with may be stated as follows: 
A unit e.m.f. is directly impressed on the terminals of a transmission 
line of length s, the distant end of which is closed by a condenser C . 
The line is supposed to be non-dissipative, its constants being in- 
ductance L and capacity C per unit length. Required the current at 
any point x(x<s) of the line. 

We write -\/L/C = k, l/\/LC = v: then by virtue of the preceding 
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analysis of transmission line propagation the indicial admittance a x 
of the infinitely long line is given by 

a s = 0, for t<x/v, 

= -r, for/^.v/i/. 
k 

The operational characteristic impedance is, of course, k = y/L/C, and 
the terminal impedances Zi and Z 2 are given by 

Zi=0, 

Z 2 = l/pC . 
Referring now to equation (244) we have: — 
X = l, Mi = l, 

_ k-l/pCo kC p-l 

M2 ~ k+i/pc ~ kCop+r 

Consequently, referring to equations (247), we have, operationally, 

v = 1 

kCop-l 



Vi =Vz = 



v 3 =v 4 = 



V:, = IV, = 



kCop+l 

kCop-l \- 
kCop+l 

kCop-l 
kCop+l 



In order to determine these functions we have therefore to solve the 
general operational equation 

kCop-l\" 



V "~\kCoP+l) 



p+ 

where V„ denotes either i'2»-i or i' 2 „. 

In order to eliminate the coefficient kC , we make use of theorem 
VIII, and write 

(p— 1\" 



MS)- 



^P+ 
In accordance with that theorem 

V„(t)=<t>n(t/kC ). 
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We therefore start with the operational equation 



MS)"- 



\p+\ 

Now the solution of this operational equation is very easy and can 
be expressed in a number of ways. We require it expressed in the 
form most easily computed. The following appears best adapted 
for our purposes. Consider the auxiliary operational equation: — 

+ ...+(- l)-2")i, 

The explicit solution is gotten by replacing 1/p" by t"/n\ and p" by 
d H /dt n , whence 

^) = (^- 2 n^+--+(- 1 )" 2 ")fr 

- i _ « ^ w (n-l) (2Q» - r- 1 V @fi! 

~ 111!" 1 " 2! 2! T " ' V ' »I ' 



But writing 

mp) 



**-[ P ) ~H(i 



it follows that 



>B V/H-l/ Jf(£+1) 

£+1 j> 1 

/> >+lH(/>+l) 

1\ £ 1 



= 1 + 



pj p+iH(p+iy 

Referring now to theorem VII, we see that 

0„(/) = (l+ j l dt}a„(t) -e-'. 

Since we have already solved for o- H (/), this determines <£„(/) and 
hence 7 W (0- The functions v 0l V\,v% ... are therefore determined. 
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Now refer back to equation (248) giving the required current in 

terms of v ,v h V2 . . . and the admittances a x (t),a« s - x (t) It 

follows at once by substitution of the preceding that 



the functions v a , v\, v?. being zero for negative values of the argument. 
This result may possibly require a little explanation. 
Consider the expression 



*jftt-rMM* 



where 1(/) denotes a function which is zero for t<t and unity for t>t . 
It is evidently identical with the admittance a x (f) provided the 
proper value is assigned to t . 

Now since l(t)—0 for t<t and unity for t^.t , the preceding may 
be written as zero for /</ , and 



m£ f{i - T)dT !ocl -'° 



which is equal to f(t — t ). 

If we set x = 0, we get the current entering the line; thus 



A {t) = Uv {t)+v x (t- |) +vt(t- f 



2s \ 



■i(»+»^-T)+^('-7) 
+»F.(.-t)+ .-..). 

This has been computed for the case where ■\/sLC = 10 and is shown 
in Fig. 2G. Referring to this figure we see that the current jumps 
at t = to the value -\/C/L = l/k, and keeps this constant value for a 
time interval 2s/v. At this instant the first reflected wave arrives and 
the current takes another jump, of 2/k. Thereafter it begins to 
decrease very slowly until time t = £s/v at which time it takes another 
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jump of 2/k. Thereafter we have a series of jumps of 2/k at time 
intervals 2s/v, the current decreasing between successive jumps. 
The smooth curve is the indicial admittance of an oscillation circuit 
consisting of an inductance sL in series with a capacity C . We see 
therefore, that the current in the line oscillates with discontinuous 




Fig. 26 — Current entering non-dissipative line terminated by capacity C„ unit E.M.F. 

applied to line 



jumps about the current in the corresponding oscillation circuit. 
Since the whole circuit contains no resistance, the oscillations never 
die away, but continue to oscillate, as shown, about the curve 



^ s Hvkc) 



which is the indicial admittance of the corresponding oscillation 
circuit. 

I shall now discuss a method of solving circuit theory problems, 
quite generally applicable to complicated networks, and particularly 
useful in dealing with transmission lines terminated in impedances. 
I have found it particularly useful in arriving at numerical solutions 
where other methods prove far more laborious. It is also of mathe- 
matical interest, as it applies another type of integral equation to the 
problems of electric circuit theory. 

Suppose that we have a network with two sets of terminals as shown 
in Fig. 27. 7 Now suppose that terminals 22 are short circuited and a 
unit e.m.f. inserted between terminals 11. Let the resultant current 
flowing between terminals 11 be denoted by Su(t)=Sn and that 

7 Regarding conventions as to signs, see the Appendix to this chapter. 
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between terminals 22 by Su(O s 5si. 5u is the driving point indicial 
admittance with respect to terminals 11 and S 2 i the transfer indicial ad- 
mittance of terminals 22 with respect to 11 under short circuit con- 
ditions. 

Similarly if terminals 11 are shortcircuited and a unit e.m.f. inserted 



Network 



Fig. 27 



between terminals 22 the current flowing between terminals 22 is 
denoted by Stt(t) =S«2 and that flowing between terminals 11 by Sn(t) 
= 5i2. If the network is passive, i.e., contains no internal source of 
energy, it follows from the reciprocal theorem that S 2 i = S u . As far 
as the two sets of terminals are concerned, the network is completely 
specified by the indicial admittances Sn,S22,S2i = Si2. 

Now let a voltage Vi{t) = V\ be inserted between terminals 11, and 
a voltage V»(J) = V 2 between terminals 22. The current flowing 
between terminals 11, denoted by I x is 

Ii(t)=j [ J^Wr)S 11 (t-r)dr+j t fWr)S 12 (t-T)dr (253) 

while the corresponding current between terminals 22 is 

T 2 (f) =jJWr)S 22 (t^r)dr+f t fWr)S 2l (t-T)dr (254) 

Now consider two networks of indicial admittances Sn, S22, Si2 = Sai 
and Tu,Tn t Tu = T21 respectively and let them be connected in tandem 
as shown in Fig. 28 to form a compound network. 

_I,3 . , 2 



?— 



S U' S 22' S I2 




T \f T V? T \2 



"I 



Fig. 28 

We require the indicial admittances of the compound network in 
terms of the indicial admittances of the component networks. 

Short circuit terminals 22 of the compound network and insert a 
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unit e.m.f. between terminals 11. Let V 3 \{t) denote the resultant 
voltage between terminals 33 measured in the direction of the arrow, 
and I 3 the current flowing between the networks. We have then the 
two following expressions for the current la. 

I 3 = S 21 (t)-j t f l V 31 (r)S 22 (l-r)dT (255) 

and 

I 3 =j t £Vn(T)Tu(t-r)dT. (256) 

Equating we get 

j t fVn{r)\S, i {t-r) + T ll {t-T)]dr = Sn_ 1 {t). (257) 

By precisely similar reasoning, if terminals 11 are short circuited 
and a unit e.m.f. inserted between terminals 22, and the correspond- 
ing voltage across terminals 33 denoted by F32, we have 8 

jJVsM [S n (t-T)+Tu(t-r)]dT = TM. (258) 

Equations (257) and (258) are integral equations of the Poisson 
type which completely determine Vn and F32 in terms of the indicial 
admittances 5 and T. We shall discuss the solution of these equations 
presently. 

If Un,V~22,U2i=Ui2 denote the indicial admittances of the com- 
pound network we have at once 

Un = S n (t)-j t f l V 31 (r)S 1 2(t-T)dT (259) 

U22 = T22{t)-j t £v 3 2{r)T2i{t-r)dr (260) 

U2l=U 1 2=J t f'V 31 (T)T 2 l(t-T)dT 

= j t fVz2{r)S l 2{t-r)dr. 



and 



(261) 



If, therefore, equations (257) and (258) are solved for V 3 i and F 32 , 
the required indicial admittances of the compound network are given 

8 V 3 -2 being opposite to V 3l in direction. 
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by (259), (260) and (261) in terms of the indicial admittances of the 
component networks. 

A simple example will now be worked out illustrating the method 
of solution just discussed. Suppose that a unit e.m.f. is impressed 
on a transmission line (infinitely long) of distributed constants R,L,C, 
through a terminal resistance R - Required the terminal line voltage 
V. 

The operational equation of this problem is gotten in the usual 
manner. The current entering the line is 



VLp+R' 

It is also obviously equal to — (1— V) : equating the two expressions, 

Ro 



and rearranging we get: — 

\Lp+R 



v= • <* 



Writing i?/2L=p and setting R = \'L/C, this becomes 

vi+m_ (262) 

l + Vl+2p/p 

This operational equation can, of course, be solved in a number of 

ways, though, as a matter of fact, its numerical solution is quite 

troublesome. This point will be returned to later: we shall first 

formulate the problem in accordance with the method just discussed. 

The indicial admittance of the line is known; it is 



G- 



r<I (pt)=A(t). 

Consequently the current entering the line is explicitly 

d <-' 



±f Q V(T)A(t-r)d, 



But the current is also equal to -=- (1— V(t))\ equating, we get 

Ro 

d r< 



V{t) = l-R j t f'v{r)A{l-r)dr. 
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Performing the indicated differentiations 

V(t) =l-RoA{o) V(t) - Rof V(r)A'(t- T )dr. 

Nowi4(o)= J^and 

A'(t)=pe-?'(I 1 ( P t)-Io(pt))JC. 

and R = \/L/C; therefore the equation becomes 

V(t) = \ + lfj(t - r) [J„ (pr) - h (pr) \e~^dr. 

As a matter of convenience we change the time scale to pt, and get 

V{t)= l 2 + \£v(t-T)[I (T)-h{T)]e-*dT (263) 

=\-\£ dTV «-4r TI °^ 

where it is understood that t is actually pt. This is the integral equa- 
tion of the problem and is in the canonical form of Poisson's integral 
equation. 

Before solving this equation numerically I shall show how a simple 
approximate solution is obtainable immediately; an advantage often 
attaching to this type of integral equation. 

The function — e~ t I (t) is equal to —1 for / = and converges 
at 

rapidly to zero. V(t) has, as we know from the operational equation, 

the initial value 1/2 and the final value 1. Neither function changes 

sign. It follows from the mean value theorem that the equation 

can be written as 

where a^l. Integrating 



and 



V{t)=\-\v{t)[e-"I {*t)-l] 



V ^= l + e-W < 264 > 
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The correct initial and final values of V{t) result for all final values 
of a<l; so that approximately 

1 



V{t) = 



l+e-'I (t)' 



This equation, while not exact, except for / = and t very large, shows 
faithfully the general character of V{t) and the way it approaches its 
final value unity. For large values of / 

whence 

V(t) = J-7=, t > 8. (264-a) 

Approximations of the foregoing type are not always possible and 
may not be of sufficient accuracy. I shall therefore give next a 
method of numerical solution which is generally applicable to integral 
equations of this type and works quite well in practice. We shall 
write the integral equation in the more general form 

u{x) =f(x) + f X u(x-y)k(y)dy (265) 

where f(x) and k(y) are known and u{x) unknown. The method 
depends on the numerical integration of the definite integral. Let 
us divide the x scale into small intervals d and for convenience write 

u{nd) =u„ 

f(nd)=f n 

k(nd) =k„. 

Now from the integral equation we have at once 

u(p)=iio=fo, 

u(d)=ui=fi+ I u(d — y)k(y)dy. 
Jo 

Now if d is taken sufficiently small 



i 



whence 



'd d 

u(d-y)k(y)dy = -^-[iiiko+Ucki], 



«i =/i+ -?r [ihko+Uoki] 



and 

1 



Ul= l^kj/2 Ul + U ° kld/2] 
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which determines Hi since u is known. Similarly 

U2=f-i+d [^Uok 2 -\-Uiki-\-h'hk ] 
which determines u*. Proceeding in the same manner 

ti*=fa+d \\iiok3-\-Uik2-\-Uiki-\-\tiiko\, etc. 

In this way we determine the value of u(x), point by point from the 
recurrence formula 

_fn-\-d \\Uok„-\-Uik n -\ + U-ik„-2+ . . -\-U„-\ki] f€jaa\ 

The result of the application of numerical integration, in accordance 
with formula (2G6), to the integral equation (263) is shown in Fig. 
(29). The dotted curve is a plot of the approximate solution as 
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Fig. 29 — Line terminal voltage unit E.M. F. i in pressed on line through resistance 

r„ = Vl/c 



given by equation (204), for a=l. We see that the voltage starts 
with the value 1/2 and slowly reaches its ultimate value, unity, 
its approach to unity, for large values of /, being in accordance with 
the formula 

V(t)= 1 —= . 

1 + 1/\/2tt/ 

The application of the foregoing method to the transmission line 
problem proceeds as follows. Let Suit), Sw{t) and Sn(t) be the 
short indicial admittances of the line. Sn(t) is the current entering 
the line (at x = 0) with unit e.m.f. directly impressed and the distant 
end short circuited. Su(t) is the current at x = s under the same 
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circumstances. Consequently from (252) 

S n (t) =a (t)+2a 2s (t)+2a is (t)+ • • ■ 

(207) 
Si2(t)=2{a s (t)+aUt)+aUt)+ ■ ■ ■<■ 

S 22 is clearly equal to Sn by symmetry. 

Now suppose that an e.m.f. E = E(t) is impressed on the line at 
.v = 0, / = 0, through a terminal impedance Z i} and the distant end 
(x = s) closed through an impedance Z 2 . We suppose these terminal 
impedances and the actual impressed e.m.f. replaced by the actual 
line voltages V\ and V 2 , impressed directly on the line at ac = and at 
x = s are 



m=j t £su{t-T)v x (i)dT 

Ii{t)—2j£SiM(t-T)Vi(r)dr 

+ifM- T)v * {r)dT - 



(268) 



(269) 



But the current at x = s is also equal to the current in the terminal 
impedance Z 2 in response to the terminal voltage F 2 : denoting by 
a 2 (t) the indicia! admittance of Z 2 it is 

Is(t)=£ t £a 2 {t-T)V 2 (r)dT. (270) 

Similarly the current entering the line at x = () is the current flowing 
in the terminal impedance Zi in response to the e.m.f. E— V\. De- 
noting by «i(/) the indicial admittance of Z x , it is 

I»(t)=i f' ai (t-r))E(r)- Wir) \dr. (271) 

(it ,/o 

Equating equation (208) and (271) and (209) and (270) we eliminate 
Io(t) and /,(/) and get 

/'[5„(/-r)+«,(/-r)l \\(j)dT- l's i2 {t-r)V.,(r)dT 
./o Jo 

= I ai(t-T)E( T )dT, (272) 

.70 

- f i S n {t-T)Vi{r)dT-\- l'[S 2 . i (t-T)-a 2 (l-T)]V 2 (T)dT = i). (273) 
Jo .70 
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These two equations are simultaneous integral equations of the Poisson 
type in Vi and Vz, which completely determine these voltages provided 
the admittances and the impressed voltages are known. They therefore 
represent the application of a new type of integral equation to the problem 
of electric circuit theory. 

The numercial solution of the general case, either by (248) or 
(272-273) is necessarily laborious when the terminal impedances are 
complicated and is only justified when the technical importance of the 
problem is considerable. I wish, however, to emphasize two points 
in this connection : the numerical solution is always entirely possible 
and, compared with other and older forms of solution, enormously 
simpler. One has only to inspect the classical forms of solution of 
problems of the type to realize the truth of this last statement. 

I shall now give two applications of equations (272-273) to specific 
problems, in one of which an approximate solution of the integral 
equation can be gotten, and in the other of which numerical integra- 
tion is applied. 

Problem I. Given a non-inductive cable of distributed constants 
C and R and length s, with unit e.m.f. applied at .v = 0, while at x = s 
the cable is closed by a condenser C a . Required the terminal voltage 
V{t) across the condenser C . 

We first write down the short-circuit indicial admittances of the 
cable; from equation (168) of a preceding section and equation (267) 
they are : — 

\—n~ { \(S 16/3 _3(tfJ | . _ . v 

= \l ,Rt ( l+2<TT+2*--+2 e .+..}. (274) 

S l2 (t)=S 21 (t) 

-»J£M-*-*+.-*+..}. (275) 

where /3=5 2 i?C/4. 
Now the current at x=s is equal to 

S 12 (t)-j t £v(r)S 2 ,(t-r)dT. 

It is also the condenser current due to the voltage V(t) ; that is 

c 4 m - 
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Equating the two expressions and integrating we get 

C a V{t) = CsMdr- f l V(r)S M (t-T)dr (276) 

,/o «/o 

which is the integral equation of the problem. In order to get an 
approximate solution without detailed computation we assume that 
the cable is long. In this case the leading terms of (274) and (275) 
are large compared with the terms following: Furthermore Sn(t) 
builds up very slowly while 5 , 22 (/) is a rapidly varying function. A 
good approximation therefore results if we take V(t) outside the 
integral sign in (276) and write 



CoV(t)= f's vl (T)<h-V(i) fs,.(T)A 



whence 



. f'SitW 

' Jo 

1+i \ H Sn(t)dt 



V{t) = ±- *? „ -. (277) 



This approximation is quite good for long cables and shows the way 
V(t) builds up quite truthfully. We see that V is initially zero, and 
builds up ultimately to unity. For large values of /, it becomes 

/ S n {t)dt 
V(t) = -,r, • (278) 

/ S ft (t)dt 

Jo 

This is the approximate formula also for the open circuit voltage, as 
may be seen by setting C o = in (277). 

In electric circuit problems, it is often sufficient, as implied above, 
to know qualitatively the behavior of an electric system without going 
through the labor of detailed computation. For this purpose the 
formulation of the problem as a Poisson Integral Equation is par- 
ticularly well adapted. A simple example will be given, which can 
be checked from the known solution. Suppose that we require the 
voltage V at point x of an infinitely long transmission line (L,R,C) 
in response to a unit e.m.f. impressed at x — 0. This is, of course, 
known from formula (211-a) : we shall here be concerned, however, 
with approximate solutions from the Poisson integral equation of 
the problem. 
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If a. v (0 denote the indicial admittance of the line at point x, then 
the current at point * is simply a x (t), which is given by formula 
(210-a). Bui if V(t) is the voltage at point x, the current is also 
given by 

jjfV(rK(/-r)rfr. 

Equating these two expansions, we get the integral equation of the 
problem 

^ f , V{r)a {t-T)dr=a s {j). 

at, Jo 

xR \~c 
Now if we write T = pt—A where p=R/2L and A = L ^-^_ i then 

ax = ylj-e-<- T + A VoVT(T+2A), T>0, 

and in terms of the relative time T, the integral equation is reducible to 

-^ f r V(T-T)e- T I (T)dT = e-^+ A U (VT(T+2A) 
(IT Jo 

while the exact formula for V is by (211-a) 



V ' Jo Vt(t+2A) 

From the integral equation it is easy to establish superior and inferior 
limits for V(T) ; it is 

Jo\J ) 

C T e-T l (T)V a {T)dT 

^ f - ~ = V b (T). 

/ e~ T I {T)dT 
Jo 

Both formulas give the correct initial and final values of V] namely 
e~ A and unity. Since V lies between V a and Vb, the mean value 
(V a -\-Vb)/2 also has correct initial and final values and should be a 
better approximation than either. The table given below shows the 
orders of approximation obtainable from the case where A =3. It 
is evident from this table that the foregoing approximate formulas 
exhibit the form of V{T) qualitatively in a quite satisfactory manner. 
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v„ 


v„ 


UV a +V„) 


V 





0.05 


0.05 


0.05 


0.05 


2 


. 25 


0.12 


0.18 


0.17 


4 


0.39 


0.19 


0.29 


0.26 


6 


0.50 


0.23 


0.36 


0.32 


8 


0.57 


0.27 


0.42 


0.37 


10 


0.64 


0.31 


0.47 


0.41 


12 


0.69 


0.34 


0.51 


0.44 


15 


0.74 


0.38 


0.56 


0.48 


18 


0.78 


0.41 


0.60 


0.52 



Problem II. Our second illustrative problem may be stated as 
follows: — A unit e.m.f. is impressed on a transmission line of length 5 
and distributed constants L,R,C. At x=s the line is closed by a resist- 
ance R in parallel with an inductance L . Required the current in 
the terminal resistance. If V{t) denotes the terminal voltage, the 
current at x = s is given by 



SM-j l f o 'v(T)S 22 (t-r)dT. 



It is also equal to the current flowing into the terminal impedance; 
that is 

i-V(t)+± f'v(T)dr. 
■Ko -LoJO 

Equating and rearranging 

[j£ +522(0) ] V{t) = ^2W- t £V(r)[£+5' 22 (/-r)]rfr. (279) 

Now the short circuit admittance S 22 and 5 i2 are given by formula 
(210-a) of a preceding chapter, and S 22 (o) =\/C/L. In order to apply 
numerical integration to (279), numerical values must be assigned to 
the constants. We take 

R„ = \/ L/C =1935 ohms, 
L o = 0.4 henry, 

|=P = 292, 
w = 1/\/LC=1.105X10 4 , 

5 = 100. 
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The results of the numerical evaluation of equation (279), with these 
values inserted, is shown in Fig. 30. The voltage is identically 
zero until *tf = 100; t = \W/v is the time of propagation of the line. 
At that instant it jumps to the value e~ pl = e ~ l00p/v and then begins 
to die away rapidly due to the draining action of the inductance. 



* 



105 110 

Values or vt 



Fig. 30 — Voltage across terminal impedance on smooth line 

The effect of secondary reflection is insignificant and therefore not 
shown. The current in the terminal resistance is V/R so that it is 
given by the same curve. 

I have reserved until the last the exposition of the expansion theorem 
solution as applied to transmission lines with terminal impedances, 
for the reason that it is the least powerful and the most restricted, 
although most closely resembling the classical form of solution. Fur- 
thermore, it does not represent the sequence of physical phenomena, 
in fact it is not a wave solution, but a solution in terms of normal 
or characteristic vibration. In practical application its usefulness is 
restricted to the non-inductive cable. 

It will be recalled that the expansion theorem solution is formu- 
lated as follows : — 

If A=l/Z(p) 

is the operational equation of the problem, then the explicit solution is 

M 

e Pkt 



*»-2Gi + 2 



pkZ'ipk) 



where p lt p 2 ... are the roots of the equation Z(p)=0. 
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Let us apply this formula to the case of a line of length s, with unit 
e.m.f. directly applied at x = s, and line short circuited at x=s. Re- 
ferring to equation (244) and putting X = ) u 1 = M2 = l we get 

l CMhTfr-.) . 1 (280) 

A sinh 7^ A(/>) 

as the operational formula of the problem. This can he written as 

^ > = (C/> + G) C -^,^ ) - JL: (281) 

r 7 smh 7^ A(/>) 

where in the general case, 

y = V(Lp+R)(Cp + G). (282) 

The values of 7 for which Z x (p) vanishes are the roots of the trans- 
cendental equation 

sinh 75=0 

excluding zero. These roots are infinite in number: Let y m be the 
m tu root; then 

y n ,=i~, m-1.2, ...00. (283) 

The corresponding values of p m are then gotten by solving (282) for 
p and writing 7 = 7™. 

The explicit solution of the operational equation (281) is then 



^»-zfe+2( c+ s) 



G \ c osh y„,(s-x) ^ 
^7* 



57m -,— cosh y,„s 

dpm 

(284) 

Av(o) ^* __, «7»i 

rfPm 

Let us apply this to the non-inductive, non-leaky cable in which 
L = G = and 7 = VRCp, so that 



s-RC 



and 



</7>» _ RC 
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Also Z x (o) = sR. We thus get 

. . . 1 , Z \"V W/.7T -To?- 1 , nnr s 

^- ( ' )= ^+rij2. ros .r"' ■ (285) 

#w=l 

This is a thoroughly practical formula for computation, owing to 
the rapid convergence of the series. In fact, for this particular line 
termination chosen, it is probably the simplest and most easily com- 
puted form of solution. These advantages depend, however, strictly 
on two facts. First, the fact that the line is taken as non-inductive 
and secondly that the terminations chosen are those of a short circuit. 
In fact, as we shall see, it is only in the case of the non-inductive 
cable that this type of solution is of any practical value. 

There is one other point which should be carefully observed in 
connection with this solution (285). This is that it is not expressed 
in terms of a series of direct and reflected waves, corresponding to the 
sequence of physical phenomena, but in terms of normal or characteristic 
vibrations. This point will be returned to later. 

Let us now attempt to apply this type of solution to the trans- 
mission line, L,R,C,G. Writing 

p ~2L r 2C 

_R L _G_ 

a 2L 2C 



We have 



whence 



v = l/VLC. 



V" 



P»- = - p±v yjyl+£ 



= -p±w^{— ) --£, m-1.2, 



dy„, 1 t , v 




i \ fmir\ 2 


V 1 ' 
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Setting G = for simplicity and substituting in (284) we get, after 
easy simplifications, 



cos 



(r-) 



^-a+^2 i/x. . "("V0r)'-g) '""• (28,i) 

\'Vt] ~* 

If we write 

|/W7T\ 2 p 2 mir 

(280) can he written as 
^*(0 = -d+ — / ) p,n (p,H?'/-.v)+sin — (nmvt+x) -. (28/) 

.S'A 5 -— W IHT ' S S ) 

Urn 

S 

This type of solution is often referred to as a wave solution and the 
component terms of the series regarded as travelling waves. As a 
matter of fact it is a solution in terms of normal or characteristic 
vibrations, each of which is to be regarded as instantaneously pro- 
duced at time / = 0. The solution in terms of true waves has been 
fully discussed in the preceding. 

Formula (287) is practically useless for computation on account of 
the slow convergence of the series (the series are only conditionally 
convergent), and cannot be interpreted to bring out the existence of 
the actual direct and reflected waves and the physical character of 
the phenomena it formulates. In fact, as stated above, this form of 
solution is useful only in connection with the non-inductive cable. 

In the cases considered above we have taken the simplest possi- 
ble terminations — these of short circuits in which case the roots of 
Z(p) are easily evaluated. If, however, the line is closed by arbi- 
trary impedances, the case is quite different, and the location of the 
roots becomes, except for simple impedances, and then only in the 
case of the non-inductive cable, practically impossible. While, there- 
fore, the expansion theorem solution can be formally written down, 
its actual numerical evaluation is a practical impossibility, except 
in a few cases. For this reason it will not be considered further here. 

The physically artificial character of the expansion solution, as 
applied to transmission lines, may be seen from the following con- 
siderations. When a wave is sent into the line, for a finite time 
equal to the time of the propagation of the line, it is independent of the 
character of the distant termination. Yet in the expansion solution 
every term involves and is dependent upon the impedance constitut- 
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ing the distant termination. Evidently, from physical considerations, 
the series of component vibrations making up the complete solution 
must therefore so combine as to annihilate the effect of the distant 
termination for a finite time. The solution is, therefore, mathe- 
matically correct but physically artificial. 

Note on Integral Equations. 

An integral equation is defined as an equation in which the unknown 
function occurs under a sign of integration; the process of determining 
the unknown function is called solving the equation. 

Integral equations are of great importance in mathematical physics 
and in recent years very considerable work has been done on them 
from the standpoint of pure analysis. 

The types of integral equations with which we are concerned in 
the present work are Laplace's Equation 



np) = r 

Jo 



s-"'f(t)dt 

fJO 

and Poisson's Equation 



*(*)=/(*)+ f X <Ky)K(x-y)dy. 
Jo 



But little work has been done on Laplace's Equation from the 
standpoint of pure analysis; its most extensive and useful applications 
appear to be in connection with the Operational Calculus. Practical 
methods of solution are extensively discussed in the text. 

We shall now briefly discuss the solution of Poisson's Equation. 

The formal series solution, which is absolutely convergent, is ob- 
tained by successive substitution. Thus suppose we write 

<t>(x)=4>o(x)+<l>i(x)+4>2(x)+ .... 

and define the terms of the series in accordance with the scheme 

0i(x) = £*<t>o(y)K(x-y)dy, 

<h(x) = fj<t>i(y)K(x-y)dy, etc., 

the resulting series satisfies the integral equation and is absolutely 
convergent. It is, however, practically useless for computation or 
interpretation. 
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A power series solution, when it exists, can be gotten by repeated 
differentiation; thus 

0'(s) */'(*) +*(«)£(*)+ \ X <t>'(x-y)K{y)dy, 

cf>'(o)=f'(o)+<t>(o)K(o) 

In this way all the derivatives at *=0 are calculable; let them be 
denoted by <j> , <t>i, 4>2 ■ • ■ Then 

This form of solution, also, is of limited practical usefulness, except 
for small values of x. 

A number of mathematicians, including Wittaker and Bateman, 
have studied the question of numerical solution and suggested other 
processes. After quite extensive study of the question, however, the 
writer is of the opinion that point-by-point numerical integration 
like that discussed in the text is, in general, the most practical, rapid 
and accurate method of numerical solution. This judgment is con- 
firmed by G. Prasad who, in a paper on the Numerical Solution of 
Integral Equations delivered before the International Mathematical 
Congress (Toronto, 1924), discusses the whole question and arrives 
at the same conclusion. 

In the text, numerical integration is carried out in accordance 
with Simpson's Rule. It is possible, of course, to employ more com- 
plicated and refined formulas for approximate quadrature. It is the 
writer's opinion that this is hardly justified in practical problems and 
that the required accuracy is more simply obtained by employing 
smaller intervals. 

Appendix to Chapter IX. Note on Conventions as to Signs in Networks 

In the network shown on page 196 the arrows indicate the direc- 
tions chosen as positive in the network itself, quite regardless of the 
presence of any e.m.fs. and currents. 

The sign attributed to a current, an e.m.f., or a voltage is positive 
if the current, e.m.f., or voltage is in the positive direction; otherwise 
the sign is negative. 

Stated more fully: 

A current at a specific point (at a specific instant of time) is posi- 
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tive if it is flowing in the positive direction; negative if flowing in the 
negative direction. 

An e.m.f. or a voltage between two points is positive if the potential 
increases in the positive direction between the two points; negative 
if the potential increases in the negative direction. (It may be noted 
that this convention makes the sign of a voltage the same as the sign 
of that e.m.f. which could be inserted between the two points without 
producing any effects in the network.) 



CHAPTER X 

Introduction to the Theory of Variable Electric Circuits 11 

In the preceding chapters it has everywhere been assumed that the 
networks are invariable: that is to say, that the constants and con- 
nections of the network do not vary or change with time. In many 
important technical problems, however, we wish to know, not merely 
what happens when an electromotive force is applied to an invariable 
network, but the effect of suddenly changing a circuit constant or of 
introducing a variable circuit element. In the present chapter we shall 
show that this type of problem can be dealt with by a simple extension 
of the methods discussed in the preceding chapters. 

The simplest and at the same time one of the most technically 
important problems of this type is the effect of sudden short circuits 
and sudden open circuits on an energized network or system. This 
type of problem will serve as an introduction to the more general 
theory. 

The Sudden Short Circuit 

Consider the network shown in Fig. 31. 




Fig. 31 

This network, which for generality is supposed to consist of two parts 

A and B, indicated schematically, is energized by an electromotive 

force .£(/) which produces a voltage V{t) between the points ab. 

The voltage V(t) is calculable by usual methods from E(t) and the 

constants and connections of the network, supposed to be specified. 

9 The material in this chapter is largely taken from a paper by the writer on 
"Theory and Calculation of Variable Electrical Systems," Phys. Rev. Feb. 1921. 
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We now suppose that, at reference time / = 0, a short circuit is sud- 
denly placed across ab; and require the effect of this short circuit on 
the distributions of currents in the network. The solution of this 
problem is based on the following proposition: 

The effect of the short circuit is precisely the same as the insertion at 
time / = of a voltage— V(t), equal and opposite to V(t), between points 
a and b. 

The resultant currents in the system for />0 are then composed 
of two components: — 

(1) The currents which would exist in the invariable network, in 
the absence of the short circuit, due to the impressed source E(t). 
These are calculable by usual methods. 

(2) The currents due to the electromotive force V(t) inserted at 
time / = 0, between the points a and b. These are also calculable by 
usual methods, since V{t) is itself known from the primary distribution 
of currents and charges. 

By the preceding analysis we have succeeded, therefore, in reducing 
the problem of a sudden short circuit, to the determination of the 
currents in an invariable network in response to a suddenly impressed 
electromotive force: that is, the problem to which the preceding 
chapters have been devoted. 

The Sudden Open Circuit 

The problem of a sudden open circuit in any part of a network 
can be dealt with in a precisely analogous manner, although the actual 
calculation of the resultant current and voltage distribution is mathe- 
matically more complicated. Consider the network shown in Fig. 32. 



ECt) 



lit) 



Fig. 32 

Here the network is supposed to be energized by an electromotive 
force E(t) which produced a current /(/) in the invariable network in 
branch ab. We require the effect of suddenly opening this branch. 
The solution of this problem depends on the following proposition. 
The effect of opening branch ab at reference time t = is the same as 
suddenly inserting at time t = 0, a voltage V(t) which produces in branch 
ab a current — /(/) equal and opposite to the current which would exist in 
the branch in the absence of the open circuit. 
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While this proposition is precisely analogous to the corresponding 
proposition in the case of a sudden short circuit, it does not explicitly 
determine the voltage V(t), which must be calculated as follows: 

Let the driving point indicial admittance of the network, as seen 
from branch ab be denoted by A ab (t) . Then, from the preceding propo- 
sition, it follows at once that V(t) is given by 



^j\(r)A ab (t-r)dT=-I(t), t>0. 



This is a Poisson integral equation in V(t), from which V(t) is calcul- 
able. With V(t) determined, the currents in any part of the network 
are calculable by usual methods, and consist of two components: — 

(1) The current distribution in the network due to the impressed 
source E(t) in the absence of the open circuit. 

(2) The current distribution due to the electromotive force V(t) 
inserted in branch ab at time / = 0. 

As in the case of the sudden short circuit, we have thus reduced the 
problem of a sudden open circuit to the determination of the current 
distribution in an invariable network in response to a suddenly im- 
pressed electromotive force. 

Variable Circuit Elements 

In the preceding cases of sudden open and short circuits it will be 
observed that the network changes discontinuously from one invariable 
state to another. A more general case, and one which includes the pre- 
ceding as limiting cases, is presented by a network which includes a 
variable circuit element : that is, a circuit element which varies, con- 
tinuously or discontinuously, with time. A network which includes 
such a variable circuit element will be called a variable network. Vari- 
able circuit elements of practical importance are the microphone trans- 
mitter, which consists of a variable resistance, varied by some source of 
energy outside the system; the condenser transmitter, which consists of 
a condenser of variable capacity; and the induction generator, in 
which the mutual inductance between primary and secondary, or 
stator and rotor, is varied by the motion of the latter. The case of a 
variable resistance will serve as an introduction to the general theory 
of such variable networks. 

Consider a network, energized by a source E(t) in branch 1, and 
containing a variable resistance element r{t) in branch n. The func- 
tional notation r(J) indicates that the resistance r varies with time. 
Let I n {t) denote the current in branch n, and assume that the network 
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is in equilibrium prior to the reference time £ = 0. The mathematical 
theory of this network depends on the following proposition : — 

The network described above can be treated as an invariable network 
by eliminating the variable resistance element r{t) and inserting an electro- 
motive force — r(t)I„(t): that is, an electromotive force equal and oppo- 
site to the potential drop across the variable resistance element. Conse- 
quently the current in the variable resistance branch is determined 
analytically by the integral equation 

!»(/) = ^£E(r)A ln (t-r)dr - ±£ r (r)I n (r)A nn (t-T)dr. (288) 

The first component is simply the current I {t) which would exist in 
the variable branch if the variable element were absent; hence, drop- 
ping the subscript n for convenience, the current in the variable 
branch is given by the integral equation 

1(0 = I {t)-~£r{r)I(r)A(t-r)dr (289) 

and the voltage across the variable element by 

v(t)=r(t)I(t). (290) 

Having determined /(/) and v(t) from this integral equation, the dis- 
tribution of currents in the network is calculable as that due to a 
source E(t) in branch 1 and a source v(t) in branch n of the invariable 
network : that is, the network with the variable resistance element 
eliminated. 

A very simple example will serve to illustrate the foregoing: — ■ 
Into a circuit of unit resistance, and inductance L = l/a, in which 
a steady current I is flowing, a resistance r is suddenly inserted at 
time / = 0: required the resultant current /(/). In this case we have: 

A{t) =indicial admittance of unvaried circuit 

= l-e~ at , 

r(t)=r, 

and the integral equation of the problem is : 



/(/) =/ -r J-jH (l-e-*y)Iit-y)dy 

= I -ra I I{t-y)e~ a *dy. 
Jo 
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If the solution is carried out as indicated by (291) below, and if the 
notation at = x is introduced, we get without difficulty 

I(t)=Io)l-r(l-e 1 (x)e-*) + r~(l-eo(x)e-*)-r 3 (l-e 3 (x)e-*)+ . . .}. 

where the function e„{x) is defined as : 

e H (x)= 1+z/l !+3C 2 /2l+x 3 /3!+ • • • +.v»-V(n-l)! 

= first n terms of the exponential series. 

For all finite values of the resistance increment r the series can be 
summed by aid of the identity 

C* 

1— e n (x)e-*= I dxe- x x"- l /(n-l)l 
Jo 

Substitution of this identity gives 

I(t)=I (l-r je-w>*dx\ 

l_|_ re -U+r).v 



= Io- 



1+r 



Equation (289) is an integral equation of the Vol terra type, which 
includes the Poisson integral equation as a special case. Its formal 
series solution is obtained as follows : — Assume a series solution of the 
form 

In(t)=I (t)-I 1 (f)+h(t)-h(t)+ . . . (291) 

and define the terms of the series by the scheme 

h{t)=~ £r{r)I {r)A{t-r)dr 



h+l(t)=-^£r(r)I k (r)A(t-r)dr. 



(292) 



Direct substitution shows that this series satisfies the integral equa- 
tion. Furthermore, it is easily shown that it is absolutely con- 
vergent. 

While this series solution is not, in general, well adapted for numer- 
ical calculations, it throws a good deal of valuable light on the ulti- 
mate character of the oscillations in the important case where E(t) 
and r(t) both vary sinusoidally with time. In this case, if the fre- 
quency of the applied e.m.f. be denoted by Fand that of the resistance 
variation by/, it is easy to show that the current I {t) in the unvaried 
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circuit is ultimately 10 a steady state current of frequency F. This 
follows from the fact that the definite integral which defines the current 
/„(/) is resolvable into the ultimate steady state current corresponding 
to an applied force of frequency F, and the accompanying transient 
oscillations which ultimately die away. The fictitious e.m.f. which 
may be regarded as producing the component current i\(/) is rj'(t)I (t); 
this is ultimately the product of the two frequencies F and /, and 
therefore resolvable into two terms of frequency F-\-f and F—f re- 
spectively. Carrying through this analysis, it is easy to show that 
each component current is ultimately a steady-state but poly-periodic 
oscillation, as indicated in the following table: 

Component Current Frequency 

h F, 

/. F+f,F-f, 

h F+2f,F,F-2f, 

h F+3f,F+f,F-f,F-3f, 

I x F+4f,F+2f,F,F-2f,F-4f. 

It is of importance to observe that the component currents involve, 
from a mathematical standpoint, multiple integrals of successively 
higher orders, the nth component I„(t) involving a multiple integral 
of the w th order with respect to I (t). Consequently the successive 
currents require longer and longer intervals of time to build up to 
their proximate steady-state values, so that the time required for 
the resultant steady-state to be reached cannot be inferred from 
the time constant of the unvaried circuit. 

From the preceding table it will be seen that the ultimate steady- 
state current is obtained by rearranging the series Io+Ii + h and is 
of the form 

+ 00 

^ A„ cos (Q+nu)t+B H sin (fi+wco)/ 

« = — oo 

where Q = 27r.Fand w=2ir/. 

It is interesting to note that this series comes within the definition 
of a Fourier series only when F = or an exact multiple of/. The 
steady-state solution is of very considerable importance and is con- 
sidered in more detail in a succeeding chapter. 

From the foregoing we deduce an outstanding distinction between 
the variable and invariable networks. In the latter the currents are 

10 It hardly seems necessary to remark that the reference time t=o is purely arbi- 
trary and that the resistance variation may start at such a time thereafter that /„(/) 
may be regarded as steady slate during the entire time interval in which we are in- 
terested. Going farther, if we confine our attention to sufficiently large values of /, 
the whole process may be treated as steady state. 
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ultimately of the same frequency as the impressed e.m.f., whereas 
in the former they are ultimately of an infinite series of frequencies. 

In the preceding example, the variable impedance element is a 
resistance r(t). If the variable element is taken as an inductance 
\{t) the voltage, corresponding to equation (290) is 

|x(/)i W . 

The case of a variable capacity element is handled as follows: 
Let 1/C = S and assume that S is variable: thus, S = S + <r(t). The 
drop across the variable condenser element is then 



v(t)=<r(t) fl(t)dt. 



Similarly a variable mutual inductance n(t) between branches m and 
n produces the voltages 

in branch m, and 

in branch n. This case may be illustrated by: 

The Induction Generator Problem 

In a sufficiently general form, this problem, which includes the 
fundamental theory of the dynamo, may be stated as follows : 

Given an invariable primary and secondary circuit with a variable 
mutual inductance Mf(t) which is an arbitrary but specified time 
function, and let the primary be energized by an e.m.f. E{t) im- 
pressed in the circuit at the reference time / = 0: required the primary 
and secondary currents. 

In operational notation the problem may be formulated by the 
equations: 

Znh-pMf(t)I 2 =E(t), 

-pMjXt)I 1 +Z 2 J 2 =0, 

in which Z n and Z 22 are the self impedances of the primary and sec- 
ondary respectively; Mj{t) is the variable mutual inductance; E(t) 
is the applied e.m.f. in the primary; and p denotes the differential 
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operator d/dt. By aid of the fundamental formula these equations 
may be written down as the following simultaneous integral equations : 

hit) =^£dyAn(J-y) (E(y)+Mj-\f(y)h(y)]) 
h(t) = M^£dyA „(/ - .v) ± [f(y)h (y)] . 

In these equations, A nit) and Aiiit) denote the indicial admittances 
of the primary and secondary circuits respectively (when il/ = 0): 
that is, the currents in these circuits in response to a unit e.m.f. (zero 
before, unity after time / = 0). We assume, of course, that they are 
known or can be determined by usual methods. 

It follows at once that the formal solution of these equations is 
the infinite series: 

I 1 (f)=X (t)+X i (f)+X i (t)+ . . . +X 2 „(/)+ . . . 

/ 2 (0 = Fi(/)+r 3 (0+nw+ . . . 

in which the successive terms of the series are defined as follows: 

Xoit) =j- t £dyA ll it-y)Eiy)=I it), 

Y 1 (t)=M~£dyA 22 (t-y) fyf(y)X (y)], 

X 2 (t)=M^£dyAnit-y)-j y [f(y) Yi(y)], 

Y*(t) =M^£dyAnit-y)~Uiy)X,iy)\, etc. 
In the light of formula 

m=§£f(y)Mt-y)dy 

the physical interpretation of the series solutions follows at once : 
Thus, Xoit) is equal to the current /„(/) flowing in the isolated primary 
in response to the applied e.m.f. £(/) ; the first component current 
Y\(t) in the secondary is equal to the current which would flow in the 
isolated secondary in response to the applied e.m.f. M{d/dt)f{t)X (t)\ 
Xoit), the second component current in the primary, is equal to the cur- 
rent in the isolated primary in response to the applied e.m.f. M{d/dt) 
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/(/) Fi(/); etc. The resultant currents are thus represented as built up 
by a to-and-fro interchange of energy between primary and secondary, 
or by a series of successive reactions. In the important case where 
the applied e.m.f. and the variation of mutual inductance are both 
sinusoidal time functions, of frequency F and / respectively, it is easy 
to show that each component current becomes ultimately equal to 
a set of periodic steady-state currents. Thus the component X is 
ultimately single periodic, of frequency F; Y\ is ultimately doubly 
periodic, of frequencies F-\-f and F—f; X 2 triply periodic, of fre- 
quencies F+2J, F and F—2f; Y 3 quadruply periodic, of frequen- 
cies F+3/, F+f,F-f,F-3f; etc. 

The Solution for the Steady-State Oscillations 

For the very important case of periodic applied forces and periodic 
variations of circuit elements we are often concerned exclusively 
with the ultimate steady-state of the system, and not at all with the 
mode in which the steady-state is approached : that is, attention is 
restricted to the periodic oscillations which the system executes after 
transient disturbances have died away. In this case, if the periodic 
variations of circuit elements are sufficiently small, the required steady- 
state is obtained in the form of a series by replacing each term of the 
complete series solution by its ultimate steady-state value; a process 
which is very simple in view of the physical significance of each term 
of the latter series. The appropriate procedure will be briefly illus- 
trated in connection with the variable resistance element. In view 
of the fact that we are concerned only with the ultimate steady- 
state oscillations, we can base the solutions on the symbolic equation 

I = Io - r Ml. (293) 

Here r(t) is the variable resistance element; I a is the current which 
would flow in the absence of the resistance variation ; and Z is a gen- 
eralized impedance of the network, as seen from the variable branch. 
Its precise significance and functional form is given below. 
We now suppose that I is given by 

L = Joe* il (real part) (294) 

= h (JoeM-hZe-™) (295) 

where the bar indicates the conjugate imaginary of the unbarred 
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symbol, so that (295) is entirely real. Correspondingly the variable 
resistance will be taken as 



K0 = -o («*"' + c- 1 '"') 



= re iul (real part) (296) 

= r cos w/. 

Here r is taken as a pure real quantity, which fixes the size of the re- 
sistance variation. No loss of generality is involved in this, since it 
merely involves referring the time scale to the zero of the resistance 
variation. 

The symbolic impedance Z, as employed in the theory of alter- 
nating currents, will depend on the frequency and is, in general, a 
complex quantity. Its value at frequency U/2ir will be denoted by 

Z(in)=Z 
while its value at frequency (Q-\-nu)/2ir will be written as 
Z(i(Q+n<a) ) =Z„. 
We now assume a series solution of (293) of the form 
J-/ +/ 1 +I 1 + ... 

where the terms of the series are defined by the symbolic equations 

r _ m T 

1\— — —y l0, 



T - ' yi) T 



(297) 
r(/) 
Z 



.Substitution shows that this series formally satisfies the equation. 

Starting with the first of (297) and substituting (295) and (29G) 
we get 

= -Ty\ 7 e , '«-^>'+7oe- ,(n + w >'+/oC , ' tn - (,,) '+JoC-' (n - w >' | , (298) 
or 

U = -\Jo\-^--+ 7 —\. (299) 

6 Zi Z._i j 

In (299) it is to be understood that the real part is alone to be retained. 
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Proceeding in a similar way with the equation 

we get 

In this way the steady-state series solution is built up term by 
term, the component currents being poly-periodic as indicated in a 
previous table. 

For sufficiently small impedance variations this method of solution 
works very well, and leads to a rapidly convergent solution. In 
other cases, however, the solution so obtained may be divergent, 
even when the complete series solution from which it is derived is 
absolutely convergent. The explanation of this lies in the fact that 
the steady-state series so obtained is the sum of the limits (as / ap- 
proaches infinity) of the terms of the complete series solution, whereas 
the actual steady-state is the limit of the sum. These are not in 
general equal; in particular the former may be and often is divergent 
when the latter is convergent. 

In view of the foregoing considerations it is of importance to de- 
velop another method of investigating the steady-state oscillations 
which avoids the difficulties in the formal series solution. The fol- 
lowing method has suggested itself to the writer and works very well 
in cases where the previous form of solution fails. It should be stated 
at the outset, however, that the absolute convergence of the solution 
to be discussed, while reasonably certain in all physically possible 
systems, has not been established by a rigorous mathematical in- 
vestigation, which appears to present very considerably difficulties. 

We start with the problem just discussed and, in view of the results 
of the formal series solution there obtained, assume a solution of the 
form: 

i N 

J = ^i IB e'W» w "+4 M r'< !! + w " (301) 

-N 
N 

= y^A m e i ^+" ,U3)t (real part). (302) 

-N 

Here the series is supposed to extend from m = +N to m=—N. 
Ultimately, however, N will be put equal to infinity. As before, the 
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bar indicates the conjugate imaginary of the unbarred symbol and 
(301) is therefore entirely real. 

If we now substitute (301) in the symbolic equation (293) we get, 
by (295) and (296), 

i-V I Amf ****** + 1*-*®+"*] =hJoe i9 -'+h Joe~ int 

2Z ^^ [ J 

Simplifying this equation and dropping the conjugate imaginaries 
gives : — 

V AneW+'^^Joe* 11 - -£- ^ A^+^+w (303) 

Finally, if we write 

Z(i(ft+>»co) )=Z m 
and 

r/Z m =h m (304) 

and equate terms of the same frequency on the two sides of the equa- 
tion, we get 

A N =-h\iA N -i 

A m = -JlmiAm-i+Am+i) (XlmKN 

(305) 

Ao = Jo-h (A- l +Ai). 

It will be observed that, by (305), starting with A N each coefficient 
is determined in terms of the coefficient of the next lower index. 
Thus : 

An=— HnAn-i 

A N -x=-h N - x (A N - 2 +A N ) 
h N -iA N - Z 



Similarly 



l — k N -\hn' 

hn -2A n -3 

1 
1— hit-iht/ 



A N - 2 =- - 

l-h N - 2 h N -i 
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Continuing this process it is easy to show that, for positive indices 
(m positive), 

A m = -h m C,„A m -i (306) 

where C„, designates the continued fraction 

1 



1 — /'»>/';>/+ 1 



1 



X — hm+lhm+Z . 



1—hsr-ihN 



The procedure for the coefficient A-„, is precisely similar. For 
convenience we write A- m =A m , Z- m — Z m , and r/Z-„, = h M . In 
this notation we get by precisely similar procedure 

A' m = -ti m C' m A' m -x (307) 

where C,„ designates the continued fraction 

1 i 



1 " III ''■ HI -\- l -I I / 7 / 



1 



l — JlN-iJlN 



We now put the index N equal to infinity and the continued fractions 
Cm and C m become infinite instead of terminating fractions. 
Collecting formulas we now have 

Am = lint (-■m-A. ,,,— i 

A ' = — k' C A ' , 

**ff| ''in*-' mf* m— 1 

and 

Ao = Jo-(h a A 1 +hU'i) 

whence 

Jo 



A = 



n, i r i " 



l-hJiid-ho'h'Q 



The coefficients are thus all determined in terms of J . 

The practical value of this method of solution will depend, of course, 
on the rate of convergence of the continued fractions. While no rigor- 
ous proof has been obtained, it is believed that they are absolutely 
convergent for all physically possible systems, but this question cer- 
tainly requires fuller investigation. Nevertheless any doubt regard- 
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ing the convergence of the solution need not prevent the use of the 
method in a great many problems where physical considerations 
furnish a safe guide. For example this method of solution, when 
applied to the problem of the induction generator, discussed above, 
leads to the usual simplified engineering theory of the induction 
generator and motor, besides exhibiting effects which the usual treat- 
ment either ignores or fails to recognize. 

Non-Linear Circuits 

In the previous examples discussed, the variations of the variable 
circuit elements are assumed to be specified time functions, which is 
the same thing as postulating that these variations are controlled by 
ignored forces which do not explicitly appear in the statement and 
equations of the problem. We distinguish another type of variable 
circuit element, where the variation is not an explicit time function 
but rather a function of the current (and its derivatives) which is 
flowing through the circuit. For example, the inductance of an iron- 
core coil varies with the current strength as a consequence of mag- 
netic saturation. The equation of a circuit which contains such a 
variable element (provided it is a single valued function) may be 
written down in operational notation 

Z/ + 0(/)=£(O, 
or 

Z/=£(/) -</>[/(/)]. (311) 

In this equation Z is, of course, to be taken as the impedance of 
the invariable part of the circuit, the indicial admittance of which 
is denoted by the usual symbol A(t). 

Equation (311) may be interpreted as the equation of the current 
/(/) in a circuit of invariable impedance Z when subjected to an 
applied e.m.f. E(t)—<f>[I(t)]; consequently, by aid of our fundamental 
formula, I(t) is given by 

/(*) = ^£a (t - y)E{y)dy - jj\ (t - y)*[I(y)])dy. 

The first integral is simply the current in the invariable circuit of 
impedance Z in response to the applied e.m.f. £(/); denoting this 
by /(/), we have 

I(t)=I (t)-j t f'A(t-y)^[I(y)]dy. 

This is a. functional integral equation, the solution of which is gotten 
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by some process of successive approximations. For example, pro- 
vided the sequence converges, /(/) is the limit as n approaches in- 
finity of the sequence 

I (t),h{t),h{t), ■■ ■ ,/»(/), 
where the successive terms of the sequence are defined by the relations: 

h{t)=Ut) -j t £A{t-y)4\I {y)\dy. 



I, l+ i(t)=Io (/) -j t £A(t-y)4>[I H (y)]dy. 

We shall not pursue the discussion of non-linear circuits further, in 
view of their mathematical complexity and their relatively specialized 
technical interest. The reader who is interested may, however, con- 
sult the writer's paper on Variable Electrical Systems, 11 for a fuller 
treatment of the subject. 

CHAPTER XI 

The Application of the Fourier Integral to Electric 
Circuit Theory 

The application of Fourier's series in electrotechnics is a common- 
place; the use of the Fourier integral, however, has largely remained 
in the hands of professional mathematicians. An outstanding dis- 
tinction between the series and the integral, from which the greater 
power of the latter may be inferred, is that the series represents only 
a periodic regularly recurrent function, whereas the integral is cap- 
able of representing a non-periodic function : in fact all types of func- 
tions, subject to certain mathematical restrictions which are usually 
satisfied in physical problems. 

Before taking up the application of the Fourier Integral to Electric 
Circuit Theory, we shall very briefly review the elementary mathe- 
matics of the series and integral; for a fuller treatment the reader is 
referred to Byerly, Fourier's Series and Spherical Harmonics. 12 

Consider a function <£(/), which in the region 0</< T is finite, single- 

11 Phys. Rev. Feb., 1921. 

12 In this chapter the Fourier Integral is approached from the view-point of its 
physical application and no completeness or rigour is claimed for the treatment. 
The mathematical theory of the Fourier integral is, of course, completely developed 
in treatises on the subject. The object of this chapter is merely to outline some of 
its applications. 
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valued and has only a finite number of discontinuities or of maxima or 
minima. In this region it can then be expressed as the Fourier series 



mere 



*(0-M«+2^» COi (T ! +jB " iln ( ? r0f (312) 

A,, = l£ r 4>(t)- cos ( 2 ^t)dt, 

(313) 
B„ = lj% (/)• sin ( 2 ™ t)dt. 

An equivalent series is 

*(*)-} Fo+^Fneos (^?/-9„) (314) 

i 
where 

F„=vzim (315) 

e„ = tan- 1 (-B„M„). 

This expansion is valid in the region 0<t<T, irrespective of the 
form of the function elsewhere. Let us, however, assume that the 
function repeats itself in the period T: that is 

*(*±*r)-*(0. ft = 1,2, 3.. .N. 

Then the expansion represents the function in the region —NT<t^NT. 
Finally if N is made infinite, the function is truly periodic and the 
Fourier series represents it for all positive and negative values of 
time. 

It follows from the foregoing that, if the Fourier series represents 
the function for all positive and negative values of time, the function 
must be periodic for all positive and negative values of time; otherwise 
the expansion is valid only over a restricted range of time. 

Now let us suppose that <f>(t) is non-periodic. For convenience, 
in connection with subsequent applications we shall suppose that it is 
zero for all finite negative values of time, that it converges to zero as 
t— >oo , and that 

/ " </»(/) dt 

exists. Such a function obviously cannot be represented by the usual 
Fourier series for all finite positive and negative values of time; it 
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can be represented, however, by the limiting form assumed by the 
series as the fundamental period T is made infinite. That is, we 
can assume that the function is periodic in an infinite fundamental 
period and this will not affect the expansion for finite positive and 
negative values of time. Proceeding in this way and putting the 
fundamental period T equal to infinity in the limit, the Fourier series 
(314) becomes an infinite integral and we get 

0(0=— / F(o>) -cos ((at-d(w))d<a (316) 

*"(«)= I [ r '</>(/) cosa^f| 2 +[ |'°°0(O sin cot df\*\ (317) 

tan 0(oj) = / 0(0 sin cotdt-r- J 0(0 cos cotdt. (318) 



where 



and 



This is the Fourier integral identity of the function 0(0 and is valid 
for all finite positive and negative values of time. 

In physical applications, particularly those to electric circuit theory, 
it is often convenient to employ exponential instead of trigonometric 
functions. The required transformation follows easily from the 
relation 

e' fl = cos 6+i sin 0, i = \/—l. 

Thus if we write 2t/T = co the Fourier series (312) is easily reduced 
to the form 

+ °° 



0(0 = N F{inwo)e hm » l (319) 



where 



as 



F(inco ) = ^f 4>(T)e-"' u ' T dT. (320) 

J Jo 

In precisely similar manner the Fourier integral (31G) can be written 
0(/) = / F{ico)-e^'dco (321) 

= *- / °°0(r)</r / e'^'-^dco. (322) 

«S7T,/o f/-W 
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Applications to Electric Circuit Theory 

Let us assume that at time /= —NT, an electromotive force E(t), 
periodic in fundamental period T, is impressed on a circuit of complex 
impedance Z(i<a), where o> denotes 2jt times the frequency. Required 
the resultant current /. 

For values of t> —NT the electromotive force (see formula (319)) 
can be expressed as the Fourier series 



where 



oo 

E(t) = y]F(intoo)e i " 03 ° l 

_ 00 

F(inu> ) = ~ / E( T )e- i "^ T dT. 
1 Jo 



The resultant current for t> —NT is therefore 

■v^a F(ino) ) ' transient oscillations 

/= 2 VrT^\ ei "" a ' + ' initiated at time 

If we are concerned with the current for values of />0, and if NT is 
made sufficiently large, the initial transients will have died away and 
the complete current for /^0, will be given by 

I=yg^l^>. (323) 

^-^Z(inco ) 

— 00 

This formula implies the periodic character of E(t) for sufficiently 
large negative values of time. If, however, E(t) is zero for negative 
values of time, we can employ the Fourier integrals (321) and (322) 
in precisely the same way and get, as the complete expression for the 
current for positive or negative values of time: — 

F^eto'dc* (324) 



= ^- / E(r)dr -yjT^d^. (325; 



The infinite integrals (324) and (325) formulate the current in the 
network, specified by the impedance function Z(ioi), in response to an 
electromotive force E(t) impressed at time / = 0; they therefore mathe- 
matically formulate, by aid of the Fourier integral identity, the funda- 
mental problem dealt with in the preceding chapters and solved by aid 
of the operational calculus. 
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No attempt will be made here to discuss the solution of the infinite 
integral (325), which is usually a problem presenting formidable 
difficulties, even to the professional mathematician. The general 
method of solution is by contour integration in the complex plane 
and the calculus of residues. By this process it has been successfully 
applied to the solution of special problems, and also to deriving some 
general forms of solution such as the expansion theorem solution. 13 
Compared, however, with the operational calculus, it has no advantages 
from the standpoint of rigour, and lacks entirely the remarkable sim- 
plicity and directness of the Heaviside method. 

In the direct solution of circuit problems, therefore, it is believed 
that the application of the Fourier integral is attended by few if 
any advantages, and presents formidable mathematical difficulties. 
On the other hand, there are certain types of problems encountered 
in circuit theory, where the Fourier integral is a powerful tool. 
These will be briefly discussed. 

The Energy Absorbed from Transient Applied Forces 

In many technical problems, the complete solution for the in- 
stantaneous current due to suddenly applied electromotive forces, 
although formally straight-forward, involves a prohibitive amount 
of labor. In yet others, the applied forces may be random and 
specified only by their mean square values. In such problems a 
great deal of useful information is furnished by the mean power and 
mean square current absorbed by the network, and to the calcula- 
tion of these quantities, the Fourier integral is ideally adapted. 
Its application depends on the following proposition, due to Ray- 
leigh (Phil. Mag., Vol. 27, 1889, p. 466), and its corollary. 

Let a function <j>(t), supposed to exist only in the epoch ()<t<T, 
be formulated as the Fourier integral 

4,(t) =— f °°|/(w) | • cos M-0(«)] da 

IT JO 

where 

/(«)={r/ <l>(t)cos(atdt~\ +T / <t>(t) sin ut dt~\ 

(T fT 

<t>(t)sin ut dt+ I (j>(t) cos iatdt. 
./o 

13 Bush, "Summary of Wagner's Proof of Heaviside's Formula." I'roc. Inst, 
of Radio Engineers. Oct., 1917. Fry. "The Solution of Circuit Problems.' (Phys. 
Rev. Aug., 1919). 
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Then 



JO 7T ,/0 

whereby the time integral is transformed into an integral with 
respect to frequency. 

A corollary of this theorem is as follows : 

If two functions <f>i(t), 4>i{t) supposed to exist only in the epoch 
Q<t<T, are formulated by the Fourier integrals 

*i(0 =— f l/i(«) I • cos l w/ -^i (»)] d "< 

7T ,70 

<t, 2 (t) =— / |/ 2 («) | • COS [«/-0 2 (w) ] rfo>, 

7T ,/0 

/ tfi(/)<M0^= — / |/i(«)l • :/2(w)l • cos{e 1 -B 3 )du. 

J0 TT Jo 



then 



The applications of these theorems to circuit theory proceeds as 
follows : — 

If an electromotive force E(t), supposed to exist only in the epoch 
0<t<T, is applied to a network of complex impedance Z(iu) = 
I Z(iu)) I e'P (u) we know from the preceding discussion of the Fourier 
integral, that the electromotive force E(t) and current /(/) are 
expressible as the Fourier integrals 



E(t) = -- / |/(w) | . cos {ut-d(u))du, 

IT Jo 



7T Jo 



(326) 



cos(oj/ — 0(oj) —0(u)))d(i). 



Z(iu)) 
It follows at once from Rayleigh's theorem that 

fw.-ifJ'W'A (327) 

Jo ttJo \Z(tw)\ 2 

Now let /„ be the current absorbed in branch n\ letz(iw) =|5(*w)|e ,tt ^ 
be the impedance of that branch and let E,,(l) be the potential drop 
across that branch. It follows at once from the corollary to Ray- 
leigh's theorem that 

W= f a °E„(t)I„(t)dt=— H ^iiL \s(i u )\coBa(u)da. (328) 

Jo 7T./0 \Z\lU)\- 
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Formulas (327) and (328) formulate the mean square current 
and mean power absorbed by the branch of the network under con- 
sideration, and enable us to calculate these quantities, even in the 
case of complicated networks, with a minimum of labor. Formula 
(327) is particularly well adapted to computation because the inte- 
grand is everywhere positive, permitting, in most problems, of easy 
numerical integration, whereas the analytical solution may be com- 
plicated. 

Formulas (327) and (328) have been applied to the theory of 
selective circuits, to the problem of interference from random dis- 
turbances, including static, and to the theory of the Schrotteffekt. 
For the details of such applications, which will not be entered into 
here, the reader is referred to the following papers. 

Transient Oscillations in Electric Wave Filters, Bell System 
Technical Journal, July, 1923. 

Selective Circuits and Static Interference, Trans. A. I. E. E., 1924. 

An Application of the Periodogram to Wireless (Burch & Bloeh- 
msma), Phil. Mag., Feb., 1925. 

The Theory of the Schrotteffekt (Fry), Journal Franklin Institute, 
Feb., 1925. 

The Building- Up of Alternating Currents 

Another application of the Fourier Integral, which may be briefly 
mentioned, is to the building-up of alternating currents in response 
to suddenly impressed sinusoidal electromotive forces. The investi- 
gation of this problem is of great importance to the communication 
engineer, since the excellence of a signal transmission system is to 
a considerable extent determined by the duration and character of 
the building-up phenomena. 

In long transmission systems the calculation of the building-up 
current as a time function is extremely complicated and laborious 
if not practically impossible. Furthermore we are usually not 
concerned with the current as an instantaneous time function, but 
rather with its envelope. The envelope of the current can be formu- 
lated and calculated by modified Fourier integrals, by the following 
process. 

Suppose that an e.m.f. E cos ut is suddenly applied, at refer- 
ence time £ = 0, to a network of transfer impedance 

Z(i.u)=\Z(io:)\e iB ^. 
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The resultant current /(/) may always be written as: 

/(O-s .y/; m J (l + p)cos(«/-.B) + <rsin («/ - 5) 



1 , E 



= 2"V(l + p) 2 + <r 2 ^ 77 ^ 7 cos(co/-5(co)-e) 



where 



= tan-^ 



1+P 



Evidently the functions p and <x, which it is our problem to deter- 
mine, must be — 1 and respectively for negative values of /, and 
approach the limits +1 and 0, respectively, as t-^oo. 

In an engineering study of the building-up process we are prin- 
cipally concerned with the envelope of the oscillations: hence with 

YVU + P^ + ct 2 . 

Our problem is therefore to determine the functions p and a and to 
examine the effect of the applied frequency co/2w and of the charac- 
teristics of the circuit, on their rate of building-up and mode of 
approach to their ultimate steady values. 

The functions p and a can be formulated as the Fourier integrals 



where 



P= - f °V«(X) +!>«(- X)] sin t\. ~ 
r .A) A 

-- f M [G,(X)-(?co(-X)]cos/X^ 

7T Jq A 

<r = - /°°[(2c(X)-r-(? w (-X)]sin/X^ 
it Jo A 

+— / M [Po>(X)-P u) (-X)]cos/X^, 

T Jo A 



&>(x) = z{iu+i\)\ ' s[n • 5 ^+^- 5 M- 
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These formulas are directly deducible from the fact that the ap- 
plied e.m.f., defined as zero for t<0 and E cos co/ for />0, can itself 
be expressed as 



-^r-coswtl 1 + — I sin /X-- . 
z L 7r Jo A J 



For important types of transmission systems, including the 
periodically loaded line, these formulas have been successfully 
dealt with and solutions of a satisfactory approximate character 
obtained. For further details, the reader is referred to a paper on 
"The Building-Up of Sinusoidal Currents in Long Periodically 
Loaded Lines" (Bell System Technical Journal, October, 1924). 

The foregoing must conclude our very brief account of the Fourier 
Integral and its applications in Electric Circuit theory; an adequate 
treatment of this subject would require a treatise in itself, and is 
beyond the scope of the present work. All that has been attempted 
is to give a very brief introduction to its significance in physical 
problems and a few of its outstanding applications in circuit theory. 
The reader who is interested in pursuing this subject further is 
referred to a paper by T. C. Fry on "The Solution of Circuit Prob- 
lems" (Phys. Rev., Aug., 1919), -which gives a rigorous discussion 
of the solution of the Fourier Integral by contour integration, to- 
gether with some general forms of solution of the circuit problem. 1 

1 It was planned to include in this paper a bibliography of the important papers 
bearing on the Heaviside operational method. This, however, has not been com- 
pleted, but plans call for its publication in the next issue. — Editor. 



